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We studied the feasibility of the measurement of Higgs pair creation at a Photon Linear Collider 
(PLC). From the sensitivity to the anomalous self-coupling of the Higgs boson, the optimum 77 
collision energy was found to be around 270 GeV for a Higgs mass of 120 GeV/c 2 . We found that 
large backgrounds such as 77 — »- W + W~ , ZZ, and bbbb, can be suppressed if correct assignment 
of tracks to parent partons is achieved and Higgs pair events can be observed with a statistical 
significance of ~ 5cr by operating the PLC for 5 years. 
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I. INTRODUCTION 

One of the most important events expected in particle 
physics in near future is unquestionably the discovery of 
the Higgs boson. The data from the ATLAS and the 
CMS experiments at the LHC and the DZero and the 
CDF experiments at the Tevatron hint at the existence 
of a light Standard-Model-like Higgs boson in the mass 
range of 115 - 130 GeV/c 2 [H-Q. If it is indeed the case, 
the discovery is expected to be declared within a year or 
so by the LHC experiments. 

In the Standard Model, the Higgs boson is responsi- 
ble for giving masses to both gauge bosons and matter 
fermions, via the gauge and Yukawa interactions, respec- 
tively, upon the spontaneous breaking of the electroweak 
symmetry. However, unlike the gauge interaction, the 
mechanism of the spontaneous symmetry breaking and 
the Yukawa interaction have been left untested. As a 
matter of fact, a Higgs doublet with its wine bottle po- 
tential, and its Yukawa coupling to each matter fermion 
in the Standard Model are mere assumptions other than 
being the minimal mechanism to generate the masses of 
gauge bosons and fermions. In other words, we know es- 
sentially nothing but something must be condensed in the 
vacuum to give the masses of gauge bosons and fermions. 
It is well known that the Standard Model cannot describe 
everything in the universe. An example is the existence 
of the dark matter which occupies about one-fourth of 
the energy density in the universe. The non-existence of 
anti-matter is another example. Since the gauge sector of 
the Standard Model is well tested, it would be natural to 
expect that some hints of physics beyond the Standard 
Model could be obtained via precise measurements of the 
Higgs boson properties. 
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The LHC experiments are likely to discover the 
Standard-Model-like Higgs boson. However, their pre- 
cision is most probably not enough to reveal details of 
the discovered particle (s) due to high background envi- 
ronments of proton-proton collisions. Thus, precise mea- 
surements of the Higgs boson properties by an electron- 
positron collider and its possible options are crucial to 
uncover its detailed properties which might go beyond 
the Standard Model. The International Linear Collider 
(ILC) has potential to study the properties of the Higgs 
boson(s) such as coupling strengths to gauge bosons and 
matter fermions including the top quark with high pre- 
cision, thereby opening up a window to physics beyond 
the Standard Model Q . 

In addition to the e + e~ collisions, high energy photon- 
photon collisions are possible at the ILC by converting 
the electron beam to a photon beam by the inverse Comp- 
ton scattering [5]. Physics and technical aspects of a 
Photon Linear Collider (PLC) as an option of the e + e~ 
Linear Collider are described, for instance, in ref. [6[. A 
schematic of the PLC is shown in FIG. [TJ The Higgs bo- 
son properties such as its two-photon decay width and 
CP properties can be studied in high energy photon- 
photon interaction and thus the PLC plays complemen- 
tary role to the e + e~ Linear Collider. It should also be 
emphasized that the Higgs boson can be singly produced 
in the s-channel process so that the required electron 
beam energy is significantly lower than that for the e + e~ 
Linear Collider. 

One of the most important observables to be mea- 
sured in the Higgs sector is its self-coupling, since it 
directly relates to the dynamics of the Higgs potential, 
i.e. the mechanism of the spontaneous symmetry break- 
ing. For example, a non-standard large deviation in the 
self-coupling can be direct evidence for strong first-order 
phase transition of the electroweak symmetry in the early 
universe Q. FIG. [2] (a) and (b) show diagrams of pro- 
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FIG. 1. A schematic of the PLC. The positron beam of the 
ILC is replaced with an electron beam. 
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FIG. 2. Diagrams including the Higgs self-coupling, (a): for 
77 — »- if if and (b): for e + e~ — >• ZHH. Higgs boson self- 
coupling occurs at red points. 



cesses which involve the self-coupling in 77 and e + e~ 
interactions. Recently, a prospect for studying the self- 
coupling at the ILC was reported. According to the 
study, the self-coupling is expected to be measured with 
precision of 57 % with an integrated luminosity of 2 ab _1 
at v^ = 500 GeV @. 

Measurements of the self- coupling at the PLC were 
discussed by several authors |9l-[llj. It has been pointed 
out that contributions to the self-coupling to the cross- 
section of FIG. [2] (a) and (b) are different and measure- 
ments in e + e~ and 77 interactions are complementary 
from physics point of view. In addition, as a Higgs boson 
pair is directly produced in the 77 interaction, required 
beam energy is lower, 190 GeV as described later, than 
that for the e + e~ interaction. This nature is important 
when considering energy update scenarios of the ILC. 

In ref. [9j], an order-of-magnitude estimation for back- 
ground processes was presented. However, the cross- 
section for the W boson pair production is 10 6 orders 
of higher than that for the Higgs boson pair production. 
The backgrounds from bbbb, bbcc, and ZZ production 
processes are also large and have the same final state as 
with the Higgs pairs for a low mass Higgs boson, which 
predominantly decays into bb. Given the situation, peo- 
ple had been skeptical about the feasibility of the detec- 
tion of the Higgs pair process at the PLC. 

In this work, we studied, for the first time, the Higgs 
boson pair creation at the PLC extensively with a pa- 
rameter set of the PLC based on an e + e _ linear collider 
optimized for the light Higgs boson of 120 GeV/c 2 and 
the same detector simulation framework as used for the 
ILC physics analysis. We report details of the analysis, 
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FIG. 3. Statistical sensitivity (Sat at) as a function of 77 col- 
lision energy. Black and red dots show the 5 k = +1 and 
5 k = — 1 cases. 



issues and prospects for the measurement of the Higgs 
boson pairs at the PLC. 



II. BEAM PARAMETERS 

In order to choose parameters for the PLC, we calcu- 
lated the statistical sensitivity S^at defined as 



\N(5k)-N sm \ 



L\tj<t(5k) - wsm\ 
y/L(rja(5K,) + ^bg^bg) 



(1) 



where, 5 k is the deviation of the self-coupling constant 
from the Standard Model. The constant of Higgs self- 
coupling A can be expressed as A = A SM (1 + 5 k), where 
A SM is the Higgs self-coupling constant in the Standard 
Model. N(5k) and A^sm are the expected number of 
events as a function of 5k and that expected from the 
Standard Model. <j(5k) and <jsm are the cross-section of 
the Higgs boson production as a function of 5k and that 
of the Standard Model, while L, 77, t?bg, and <tbg are 
the integrated luminosity, the detection efficiency for the 
signal, the detection efficiency for backgrounds, and the 
cross-section of background processes, respectively. For 
r] = 1 and ?7bg = 0, S'stat is written as 
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FIG. [3] plots 5 sta t as a function of the center of mass en- 



ergy of the 77 collision (denoted 



hereafter) for the 



Higgs boson mass of 120 GeV/c 2 with the 77 integrated 
luminosity of 1000 fb _1 . The cross-section of the signal 
was calculated according to the formula described in ref. 
[loj for the case of 5k = +1 and 5k = — 1 as indicated in 
FIG. [U From this result, we found the optimum energy 
to be w 270 GeV. 

The parameters for the electron and the laser beams 
are summarized in TABLE [U It was designed to maximize 
77 luminosity at <sjs^ ^ 270 GeV based on the TESLA 
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TABLE I. The parameters of electron and laser beams based 
on TESLA optimistic parameters. The polarization of elec- 
tron beam was assumed to be 100 %. 



parameter 


unit 




electron beam energy 


E e [GeV] 


190 


# of electrons / bunch 


N x 10 10 


2 


longitudinal beam size 


a z [mm] 


0.35 


transverse emittance 


~f£ x / y [10" 6 m-rad] 


2.5/0.03 


P function @ IP 


f3 x/y [mm] 


1.5/0.3 


transverse beam size 


a x/y [nm] 


100/5.7 


laser wavelength 


Xl [nm] 


1054 


laser pulse energy 


[J] 


10 


x — 4uE e /m1 




3.42 




FIG. 4. Luminosity distribution generated by CAIN. Input 
parameters are shown in TABLE [I] 

optimistic parameters [12]. The wavelength of the laser 
was chosen to be 1054 nm, which is a typical wavelength 
for solid state lasers. The electron beam energy was cho- 
sen to maximize 77 luminosity around 270 GeV, while 
keeping the electron beam emittance and the f3 functions 
at the interaction point the same as the TESLA parame- 
ters. The luminosity distribution was simulated by CAIN 
[l3| , as shown in FIG. 2J The 77 luminosity in the high 
energy region {y/s^j > 0.8 ^/s^ x ) was calculated to be 
1.2 x 10 34 cm- 2 s- 1 . 



III. SIGNAL AND BACKGROUNDS 

FIG. [5] shows the cross-sections for various processes 
of 77 and e + e _ collisions as a function of the center of 
mass energy. The figure indicates that the 77 — >> WW 
and 77 — » ZZ processes will be the main backgrounds 
at y/s 77 = 270 GeV because the total cross-sections are 
about 90 pb and 60 fb, respectively, far exceeding that of 
77 —> HH, which was calculated to be 0.19 fb. It should 
be noted that ^/s 11 = 270 GeV is below the threshold 
of the 77 — >> it process so that it is not necessary to be 
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FIG. 5. The cross-sections of various Standard Model pro- 
cesses as a function of collision energy. Solid lines show the 
77 collision case. The red line shows the optimum energy, 270 
GeV. 



TABLE II. Branching ratios of the Standard Model Higgs 
boson with a mass of 120 GeV/c 2 . 



decay mode 


branching ratio 


H -+ bb 


0.68 


H WW" 


0.13 


H^gg 


0.071 


H — >■ TT 


0.069 


H —> cc 


0.030 


H ZZ* 


0.015 


H — ► 77 


0.0022 


H -> 7Z 


0.0011 


H ss 


0.00051 


H — > 


0.00024 



considered as a background source. 

TABLE [Til shows the branching ratios of the Standard 
Model Higgs boson with a mass of 120 GeV/c 2 0. Since 
the main decay mode of the 120 GeV/c 2 Higgs boson is 
H — » bb, we concentrated on the case where both Higgs 
bosons decay into bb in this analysis. This implies that 
the 77 — >• bbbb process must also be considered as a pos- 
sible background process. 

The numbers of events expected for the signal and the 
backgrounds were calculated from the 77 cross-sections 
by convoluting them with the luminosity distribution, as 

/dL 
^{W^-^—dW^. (3) 

We used the formula in refs. fToL fTT| for the calculation of 
77 -»• HH, HELAS for 77 -)• WW, gamgamZZ-code 



TABLE III. The detector parameters. p,pr and E are mea- 
sured in units GeV. The angle is measured from the beam 
axis. 



Detector 




Resolution 


Vertex detector 


a b = 7.0 6 


3 (20.0/psin 3/2 <9) fim 


Drift chamber 




i.i x nrVr© o.i% 


ECAL 


(te/E 


= lb%/VE@l% 


HCAL 


cje/E 


= 40%/v^E 2% 



PS [13 and HELAS for 77 ZZ, and GRACE [3 for 
77 — >> 6666. The numerical integration and subsequent 
event generation was performed by BASES/SPRING 
[l9j . With this calculation, we expect 16 events/year for 
77 — >• 1.462 x 10 7 events/year for 77 — >• WW, and 

1.187 x 10 4 events/year for 77 — >• ZZ. For 77 — >> 6666, 
5.194 x 10 4 events/year is estimated for events with 66 
mass grater than 15 GeV/c 2 . 



IV. SIMULATION AND ANALYSIS 

JLC Study Framework (JSF) [20|, [2l| was used as our 
simulation framework in this study. Pythia6.4 [22| was 
used for parton shower evolution and subsequent hadron- 
aization. For the detector simulation, a fast simulator, 
QuickSim [2lj |. was used instead of a full detector simu- 
lation in order to process the huge background samples. 

QuickSim is, however, fairly detailed and realistic: it 
smears track parameters with their correlations, vertex 
detector hits according to given resolution and multiple 
scattering. It simulates calorimeter signals to individual 
cells so as to take property into account their possible 
overlapping. The calorimeter signals are then clustered 
and matched to charged tracks, if any, to form particle- 
flow-like objects to archive best attainable jet energy res- 
olution (see ref. [23| for more details). 

The detector parameters are summarized in TABLE 
IIIII In the simulation, we assumed a dead cone of a half 
angle of 7.6° in the forward/backward region of the de- 
tector to house the laser optics, the beam pipes and the 
masking system [24j . 

We generated 5 x 10 4 Monte- Carlo events for 77 — >> 
HH, 7.5 x 10 7 for 77 WW, 1 x 10 6 for 77 ZZ, and 
1 x 10 6 for 77 — >• 6666, respectively, which are statistically 
sufficient to assess the feasibility of 77 — >> HH measure- 
ment against the large number of background events. 



IV. 1. Event Selection 

First, we applied the forced 4-jet clustering to each 
event in which the clustering algorithm was applied to 
each event by changing the clustering parameter until 
the event is categorized as a 4-jet event. We used the 
JADE clustering [25[ as the clustering algorithm. 




FIG. 6. The concept of "nsig" method. A 6-hadron is gen- 
erated at "Interaction Point" and decays at the circle in- 
dicated as "Decay of 6-hadron." Arrows represent particle 
tracks. Dotted lines are some of the extrapolated particle 
tracks towards the interaction point. 



Using four- momenta of reconstructed jets, xf s = 
H, W, Z, 66) were calculated for possible jet combinations 

as 



mm 



(Mi - Mi) 2 (M 2 - Mi) 2 



J 2ji 



>2ji 



(4) 



where, M\ and M2 are invariant masses of two jets. Mi 
(i = H, W, Z, 66) are the masses of the Higgs boson, the 
W boson, the Z boson, and the invariant mass of 66 (10 
GeV/c 2 ), respectively, cF 2 ji = H, W, Z, 66) are their 
corresponding mass resolutions, and are chosen to be 
&2jH = 8 GeV and a 2jW = (?2jZ = °2jbb = 6 GeV, 
respectively. The "min[ ]" stands for the operation to 
choose the minimum out of all the jet combinations. 

In order to discriminate 6-quarks, we used the "nsig" 
method for the 6-tagging in this study. FIG. [6] illustrates 
the concept of the "nsig" method. For each track in a 
reconstructed jet, N s \ z = L/cjl was calculated, where 
L is the distance of closest approach to the interaction 
point of the track in the plane perpendicular to the beam 
and <jl is its resolution. Then N fl v (a), the number of 
tracks which have N s [ g > a, is calculated for each jet as 
a function of a. 

Before optimizing event selection criteria, we applied 
the pre-selection to reduce the number of background 
events to a level applicable to the Neural Network anal- 
ysis. The criteria for the pre-selection are 

• JVjet(JVofFv(3.0) > 1) >3, 

• JV jet (iW3.0) > 2 ) > 2 > 

• fcj > 0.05, 

• |cos6> 2 j| < 0.99, 

where, N^ et (N ^ v (b) > c) is the number of jets for which 
^Voffv(6) is greater than or equal to c. faj is the speed of 
reconstructed 2- jet which has the least xf an d 2 j is the 
angle of reconstructed 2-jet system with respect to the 
beam axis. 

After the pre-selection, we applied the Neural Network 
analysis to optimize the selection criteria. It is a three- 
layer network with a single output. JETNET [26| was 
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FIG. 7. Typical distributions of Neural Network input variables, (a): Xz after the pre-selection. (b): cos# x 2 after the W 

filter, (c): Xh after the bb filter. Black, blue, green, and red histograms show the 77 — > HH (signal), 77 — > WW, 77 — > ZZ, 
and 77 — y bbbb events, respectively. 



used to train the Neural Network system which employed 
the back propagation for the weight optimization. 

For 77 — >• WW events, inputs to the Neural Network 
are x%i x|, the visible energy, N^ t (N oSv (S.5) > 1), 
A/j e t (7V ff v (3.5) > 2), the longitudinal momentum, the 
transverse momentum, the number of tracks, and Y cut of 
jet clustering. 29958 signal events and 83777 background 
events were used for Neural Network training with the 
number of intermediate layers of 18. FIG. (a) shows 
the typical distribution of x%- 

Neural Network inputs for 77 — >• bbbb analysis are 
Xh, cos ^x|f' cos ^x^' ^ ne y i s it>le energy, the num- 
ber of tracks, Y cut of jet clusteri ng, thrust [23, pp. 284], 
sphericity [27), pp .281], Y value [23, pp .282], cosOj, and 
the largest | cos#j| of the event, where @ x 2 H {xl h ) anc ^ ®o are 
the angle of H(bb) system and of each jet, with respect 
to the beam axis. FIG. (b) shows the distribution of 
cos# v 2 after the W filter. 7756 and 1409 events for the 
signal and background, respectively, were used with the 
number of intermediate layers of 34. 

For the 77 — » ZZ events, we used Xh> Xwi Xz, 
the visible energy, the number of tracks, the longi- 
tudinal momentum, the energies of the 2-jet systems, 
JVjet(JVofFv(3.5) > 1), and 7V je t(^offv(3.5) > 2) as Neural 
Network inputs with 4536 signal and 1189 background 
events for the training with 20 intermediate layers. FIG. 
0(c) shows the typical distribution of x%- 

The Neural Network was trained to maximize statisti- 
cal significance E defined as 

E = ^signal (K) 

where, N s [ gna \ and Nbg are the numbers of remaining 
signal and background events, respectively. To reduce 
possible systematic effects from the training of the Neu- 
ral Network analysis, the performance of the Neural Net- 
work was evaluated by applying the results of the train- 
ing (weight files) to events generated separately from the 
training samples. In order to reduce the effect of the 
statistics of the event samples, we prepared the same 
number of events for the test sample for each training 



sample. TABLE IIVI shows the summary of event selec- 
tion with JADE clustering. From TABLE HV1 the statis- 
tical significance with the JADE clustering Ejade was 
calculated to be 

SjADE = 0.922^-045^ (6) 



IV. 2. Event Selection with an Ideal Clustering 

The result in the previous section indicated that it 
is necessary to improve the performance of event selec- 
tion. In order to evaluate the effect of the jet clustering, 
we applied an "ideal jet clustering" to the 77 — )• HH, 
77 — » WW, and 77 — » ZZ events where each track is 
assigned to its parent (H, W, or Z) by color information 
obtained from the event generators. The "ideal jet clus- 
tering" was not applied to the 77 — >• bbbb events since 
the color singlet combinations were non-trivial for this 
process. Input variables to the Neural Network and the 
number of intermediate layers are the same as the JADE 
clustering case. As with the previous analysis, the pre- 
selection was applied and events survived the selection 
cuts were used for the Neural Network analysis. The 
number of signal/background events used for the Neu- 
ral Network training were, 29152/57058 for 77 — >> WW, 
24305/6349 for 77 bbbb, and 22823/291 for 77 ZZ, 
respectively. FIG. [8] (a), (b), and (c) show the typical 
distributions of after the pre-selection, cos# % 2_ after 

the W filter, and Xh a fter the bb filter, respectively. We 
again applied the results of the Neural Network training 
to the event samples which are statistically independent 
of the training samples. TABLE |V1 shows the summary 
of the event selection with the ideal jet clustering. From 
TABLE the significance Sideai was calculated to be 

S id eal= 4.87 ±0.13(7. (7) 

This result indicates that 77 —> HH would be observed 
at ~ 5<j significance level with the integrated luminosity 
that corresponds to 5-year operation of the PLC, if the 
jet clustering performed perfectly. 
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TABLE IV. Cut statistics with JADE clustering. The numbers in the table are the expected numbers of surviving events 
expected in 5 years. The error on each number is from statistics of the Monte- Carlo study. 





77 — > HH 


77 WW 


77 ZZ 


77 — y bbbb 


expected events 


80 


7.31 x 10 7 


59350 


259700 


pre-selection 


47.72±0.28 


81312±282 


5172±18 


80002±144 


applying W filter 


12.27±0.14 


24.4±4.9 


231.6±3.7 


378.1±9.9 


applying bb filter 


5.867±0.097 


i.yo_ . 6 i 


59.3±1.9 


13.2±1.9 


applying Z filter 


3.766±0.078 


U -0 


5.40±0.57 


7.5±1.4 



i/) 30000 



600r 




% 50 100 150 200 250 300 350 400 450 500 

(a) %\ 



o 
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FIG. 8. Typical distributions of input variables in the case of the ideal jet clustering, (a): Xz after the pre-selection. (b): 
cos# x 2 after the W filter, (c): Xh after the bb filter. Black, blue, green, and red histograms show the 77 — »• HH, 77 — >> WW, 

77 — » ZZ, and 77 — >• fc»5fc>6 events, respectively. 



TABLE V. The similar table to TABLE El but with the ideal jet clustering. 





77 HH 


77 -> WW 


77 — >• ZZ 


77 — > bbbb 


expected events 


80 


7.31 x 10 7 


59350 


259700 


pre-selection 


46.64±0.27 


55836±233 


4172±16 


77778±142 


applying W filter 


40.13±0.25 


7.8±2.8 


46.3±1.7 


1826±22 


applying bb filter 


36.03±0.24 


7.8±2.8 


18.5±1.0 


7.8±1.4 


applying Z filter 


34.68±0.24 


4.9±2.2 


5.22±0.56 


6.0±1.2 



V. SUMMARY 

We studied the feasibility of the measurement of Higgs 
pair creation at the PLC, which is a possible option of 
the ILC. The optimum center of mass energy of the 77 
collision was found to be around 270 GeV for the Higgs 
boson with a mass of 120 GeV/c 2 . 

We found that the 77 — >> HH process can be observed 
with a statistical significance of about 5 a for the inte- 
grated luminosity corresponding to 5 years of the PLC 
running against the background process which has 10 6 
times larger production cross-section (77 — » WW) than 
the signal and other backgrounds which have the same fi- 
nal state (77 ZZ, and 77 — >• bbbb), if each track could 
be successfully assigned to parent particles (or partons). 

Our analysis showed, for the light Higgs boson, im- 
provement of the jet clustering technique is crucial to 
discriminate the backgrounds by invariant mass informa- 
tion rather than to improve the 6-quark tagging efficiency. 
This fact is reasonable, because the WW background 



turned out to be suppressed by a simple 6-quark tagging 
scheme since the W bosons do not decay into 6-quark 
pairs, while the ZZ and bbbb backgrounds can only be 
suppressed by their mass differences. 

For further improvements, vertex information from b- 
tagging analysis must be taken into account in jet clus- 
tering, thus they should be coherently developed. Efforts 
in this direction are on-going as a part of the ILC physics 
study [28|, [29j| and significant improvement could be ex- 
pected in near future. 

This analysis shows possibility to measure the Higgs 
boson self-coupling at a lower beam energy than that 
of the e + e~ mode and is useful in considering energy 
upgrade scenarios of the ILC. 
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